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Abstract 



We present a complete next-to-leading order QCD calculation for single-inclusive large- 
Pt hadron production in longitudinally polarized lepton-nucleon collisions, consistently 
including "direct" and "resolved" photon contributions. This process could be studied 
experimentally at a future polarized lepton-proton collider like eRHIC at BNL. We exam- 
ine the sensitivity of such measurements to the so far completely unknown parton content 
of circularly polarized photons. 



I Introduction 



The successful start of the RHIC spin program marks a new era in spin physics. Very inelastic pp 
collisions at high energies open up unequaled possibilities to thoroughly address many interest- 
ing questions, most importantly perhaps those concerning the still unknown role of polarized 
gluons pp. QCD analyses of upcoming high-py prompt photon, jet, heavy flavor, inclusive- 
hadron, and boson production data will greatly deepen and enrich the understanding of 
the nucleon and of QCD that we had gained previously from polarized deeply-inelastic scatter- 
ing (DIS) of leptons off nucleons. 

Investigations of interactions between polarized leptons and nucleons will, however, continue 
to play a vital role in spin physics. Besides the many observables of DIS, photoproduction 
processes are particularly interesting. They give access to the parton content of circularly 
polarized photons, due to the presence of contributions to the cross section where the photon 
"resolves" into its parton content, for instance, by fluctuating into a vector meson of the same 
quantum numbers, prior to the hard QCD interaction. Such contributions compete with the 
"direct" part, for which the photon simply interacts as an elementary pointlike particle. 

It is well-known [2^ that at lower (fixed-target) energies, the direct component tends to 
dominate over the resolved one. Only at collider energies do resolved contributions become 
important. Extensive studies in the unpolarized case at the HERA and LEP colliders have 
firmly established the concept of photonic parton densities P] and shown that resolved photon 
processes contribute a sizable fraction of photoproduction events in certain kinematic regions. 
The fact that e~^e~ and ep data can be universally described by the same set of photonic parton 
densities marks an important test and success of QCD. 

It has not yet been possible to observe resolved contributions for polarized photons, due to 
the fact that only fixed-target lepton-nucleon experiments have been performed so far in the po- 
larized case. In this context, an exciting possibility would be to have a polarized lepton-proton 
collider such as the planned eRHIC project at BNL which is currently under discussion. 

Motivated by this prospect, we present in this paper a study for photoproduction of single- 
inclusive high-pT hadrons at a future lepton-proton collider with longitudinally polarized beams, 
and investigate the sensitivity of this reaction to the contribution from resolved photons. We 
had made first exploratory studies of this, and related, reactions in Our present study goes 
beyond the results of |3] in that we now consistently include both the direct and the resolved 
photon processes in the next-to-leading order (NLO) of QCD. This is crucial for reliable quan- 
titative analyses in the future as theoretical uncertainties are expected to be under much better 
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control beyond the lowest order (LO) approximation of QCD. For the resolved contribution, 
the NLO corrections to polarized QCD hard-scattering have become available only recently in 
[HI El- For the polarized direct component, the corresponding cross sections were already cal- 
culated some time ago in jH]. We have rederived the expressions and agree with the published 
results. In addition to extending our study to NLO, we also tailor our phenomenological results 
to the eRHIC situation. 

In the next section we sketch the technical framework for the description of polarized pho- 
toproduction reactions. Section III is devoted to numerical studies for the eRHIC. We conclude 
in Section IV. 



II Technical Framework 



We consider the spin-dependent photoproduction cross section for the reaction Ip —>■ I'HX, 
where the hadron H is at high transverse momentum px, ensuring a large momentum transfer. 
We may then write the differential cross section in a factorized form: 
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This generic form applies to both the direct and the resolved cases. The subscripts and 
"H — " denote the helicities of the colliding lepton and proton, /i,, the renormalization scale, and 
S = {Pi + PpY the available c.m.s. energy squared. Neglecting all masses, V and W in Eq. ((H) 
can be expressed in terms of the c.m.s. pseudorapidity t] and transverse momentum p^ of the 
observed hadron: 



V =1- 



Pt 



and W 



Pt 



/S SV{1 - V) 

We count positive rapidity in the proton's forward direction. 



(2) 



The sum in Eq. is over all partonic channels a + b ^ c + X contributing to the single- 
inclusive cross sections, with (iAa^^^^-^ and dAa^^fj_^^j^ the associated polarized LO and NLO 
partonic cross sections, respectively. For the direct contribution to the cross section, depicted 
in Fig. n (a), we have a = 7, while for the resolved ones [see Fig. ^ (b)] a stands for the parton 
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Figure 1: Generic direct (a) and resolved (b) photon contributions to the process Ip I'HX . 

emerging from the photon. For convenience, we introduce effective polarized parton densities, 
Afj^{xi, fif), in the lepton, with Xi the fraction of the lepton momentum carried by parton a 
and /i/ the factorization scale. The A/^(?/, /x/) are defined via the convolution 



J xi y \ y J 



(3) 



with 
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y 



ln M^7^' +2mivf-^-l^)l (4) 

being the spin-dependent Weizsacker- Williams "equivalent photon" spectrum including its non- 
logarithmic contributions^ 0. Here, mi is the lepton mass and Qmax the allowed upper limit 
on the radiated photon's virtuality, to be fixed by the experimental conditions. For the direct 
case, one simply replaces 

A/J(x„/x/) 5(1 -X,), (5) 

while in the resolved contribution the yet unmeasured parton distributions of circularly polar- 
ized photons occur that are the focus of this paper. Likewise, A/^(xp, /i/) denotes the density of 
parton b in the polarized proton, and D^{z, fi'j-) is the usual unpolarized fragmentation function 
for parton c going into the observed hadron H, at scale /ij. The relations between parton-level 
and hadron-level variables are 

s.i,.^,.f....,S, (6) 

VWZ Z[l — V) 

As we mentioned earlier, it is a main feature of our paper that for the first time we include 
both the direct and the resolved contributions to the polarized cross section at NLO, as indicated 



^The latter turn out to be numerically unimportant in our study, but we keep them for completeness. 
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by the first-order corrections (iAo"^),^^-^ in Eq. IQ. The actual calculation of these corrections 
is generally a formidable task. Since later on we will mainly consider spin asymmetries, we also 
need the corresponding spin-averaged partonic cross sections. Let us go through the various 
ingredients that we need in a little more detail. 

The resolved contribution is technically equivalent to a hadroproduction cross section such 
as pp HX. The corresponding d/S.(T^^^^^-^ may therefore be immediately adopted from our 
recent calculation jB] for this case. In that reference, we also re-derived the (previously known 
[TU] ) expressions for the unpolarized case. There are in total 16 different partonic channels. 

The NLO contributions to the direct-photon component have been calculated in ^T] ^] 
and !B| for the unpolarized and polarized cases, respectively. Having the techniques available 
that we used in ^6 , it was a straightforward, though tedious, work to check the old results 
from the scratch. For the polarized case we fully agree with at an analytical level. In the 
unpolarized case, we were able to identify some minor mistakes (beyond those noted in jSJ) in 
the computer code of [12]; apart from those we agree with that calculation as well. 



The technical details of the NLO calculation for the direct part may be found in jS]. We 
will not repeat these here, but will only recall an important aspect of photoproduction cross 
sections beyond the LO: writing the cross section as the sum of the direct and the resolved 
photon contributions, 

dAcr = (iAcTdir + C^Adres , (7) 

we emphasize that neither (iAcTdir nor d/^a^cs are individually physical, i.e., experimentally 
measurable. This can be easily seen by noticing that the 2 — * 3 contributions to the direct cross 
section have singular configurations arising from a collinear splitting of the incoming photon 
into a quark-antiquark pair. These contributions are absorbed, according to the factorization 
theorem, into the so-called "pointlike" part of the NLO photonic parton densities and thus into 
the resolved component. The subtraction of the singular part is unique only up to finite pieces, 
and the ensuing ambiguity (i.e., the choice of "factorization scheme") cancels only when the 
direct and resolved parts are combined. 

As an example, when using dimensional regularization {d = 4— 2e), the appropriate "counter 
cross section" that cancels the ■j ^ qq collinear singularity in the subprocess 'jq — * gX is 
schematically given by 



dAa. 



i-counter 



where dAaqq^gg is the polarized Born cross section for the reaction qq gg in d dimensions. 
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and (cq denotes the quark charge) 



-^ + 7e- ln47r j 3eJ [2x - 1] | ) + Ahg^, 



x] 



(9) 



The freedom in the choice of the factorization prescription is reflected by the arbitrary finite 
piece Ahqy{x) which may be subtracted alongside the singular 1/e contribution. For the MS 
convention, Ahq^{x) = 0. An alternative factorization scheme (DIS^) was proposed in the 
unpolarized case JH] to facilitate the analysis of data on the DIS photon structure function 
F2 . Here the photonic Wilson coefficient 6*2,7 the NLO expression for is absorbed into 
the definition of the quark distributions of the photon. A similar scheme can be set up in the 
polarized case as well |T3] by subtracting the corresponding polarized coefficient function ACi,^ 



in the structure function gj. The photonic parton densities in the MS and DIS-y schemes are 
then related in the following way |14j : 



(10) 



where 
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The relation between the partonic cross sections for the direct contribution reads accordingly: 



dAcr 



DISt, 



(12) 



with 6Af2 given above and the symbol ® denoting a standard convolution. The dAoab^cX 
are the appropriate LO partonic cross sections. For the resolved cross section all partonic 
cross sections remain unaffected by this transformation. It is straightforward to verify that the 
scheme transformation indeed leaves the physical cross section in Eq. ((7j) invariant. 

Our explicit analytical expressions for the polarized NLO subprocess cross sections for the 
direct and resolved cases are too lengthy to be given here, but can be found in our computer 
code which is available upon request. 



Ill Numerical Results 

Let us now turn to a first numerical application of our analytical results. Instead of presenting 
a fuU-fiedged phenomenological study of single-inclusive hadron production in polarized Ip col- 
lisions, we only focus on the most interesting questions: the importance of the NLO corrections. 



5 



the dependence on unphysical scales, and the sensitivity to the yet unmeasured parton densities 
of circularly polarized photons. 

For all our calculations we choose a c.m.s. energy of = 100 GeV as planned for polar- 
ized electron-proton collisions at eRHIC. We limit ourselves to the case of neutral pions (vr*^); 
extension to other hadrons would be straightforward. We will always perform the NLO (LO) 
calculations using NLO (LO) parton distribution functions, fragmentation functions, and the 
two-loop (one-loop) expression for as- The value for as is taken according to our choice of 
proton parton distributions, which is CTEQ5 [T5] for the unpolarized case and GRSV 16 for 
the polarized one. 

In the equivalent photon spectrum in Eq. we use similar parameters as the HI and ZEUS 
experiments at HERA. We employ Q^ax — ^ GeV^ and take the photon momentum fraction to 
be within 0.2 < ?/ < 0.85. 

The fragmentation functions for neutral pions, D^, have been determined with some ac- 
curacy from analyses of data for pion production in e^e~ collisions IH?] . Even though 
our knowledge of the certainly still needs improvement, we assume for this study that 
the fragmentation functions will be well known by the time experiments at a future polarized 
lepton-proton collider are carried out. For our analysis we take the LO and NLO sets of Ref. J7] 
which also yield good agreement with recent high-p^ data for pp 7r°X at = 200 GeV by 
the PHENIX collaboration dH]. 

Since the parton distributions A/'*' of circularly polarized photons are completely unmea- 
sured so far, we have to invoke some model for them. Two extreme scenarios were considered 
in 120112], with "maximal" [A/'''(x, /iq) = P{x,ij,q)] and "minimal" [A/'^(x,yUo) = 0] satura- 
tion of the "positivity constraint" |A/'^(x, /io)| < P{,x,no). /iq denotes the scale where the 
boundary conditions [201 for the evolution are specified. For the spin-averaged densities the 
phenomenologically successful GRV photon distributions |2I] were used. We employ these two 
extreme sets for the in Eq. 0, along with appropriate LO [20] and NLO [111122] boundary 
conditions and evolution equations. Our default choice for the numerical studies will be the 
"maximal" set. 

As we already pointed out in previous LO studies for single-inclusive measurements 
rapidity differential cross sections are particularly suited for extracting the resolved contribu- 
tions and hence the polarized photonic parton distributions. This can be easily understood 
by looking at the momentum fractions of the photon and proton distributions as a function 
of rapidity. As we count positive rapidity in the proton's forward direction, large 
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Figure 2: Unpolarized and polarized 7r° photoproduction cross sections in NLO (solid lines) and LO 
(dashed lines) at ^/S = 100 GeV. The polarized results are enlarged by a factor of 10. The lower 
panel shows the ratios of the NLO and LO results in each case. 



probed at large negative rapidities. Here one expects no difference between the two photon 
scenarios since the cross section is dominated by the direct photon contribution sitting strictly 
at = 1. In addition, in the region — > 1 the Af"^ are dominated by the "pointlike" part 
which only depends on the starting scale /io but not on the details of the hadronic input. On 
the other hand, for large and positive rapidities small values are probed, and the resolved 
contribution is expected to dominate more and more. 

Figure 121 shows our results for the unpolarized and polarized photoproduction cross sections 
for neutral pions at LO and NLO at = 100 GeV as functions of the laboratory-frame 
pseudorapidity r/iab- For the asymmetric eRHIC kinematics {Ep = 250 GeV, Ee = 10 GeV), r/iab 
is related to the c.m.s. frame rapidity via 

^lab = ^ + ^ ^ ■ (13) 
We have integrated over the transverse momentum > 4 GeV of the produced pion. All scales 
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Figure 3: Scale dependence of the polarized cross section for 7r° production at LO and NLO in the 
range pr/2 < yUr = yU/ = A*/ < '^Vt- The solid and dashed lines correspond to the choice where all 
scales are set to p^- 



have been set equal: /i,. = /i/ = /^j = "Pt- 

The lower part of the figure displays the so-called "i^-factor" 

K^'^^^f^ (14) 

One can see that in the unpolarized case the corrections are roughly constant, of size 60% — 80% 
over most of the ?7iab region considered. In the polarized case, we find generally much smaller 
corrections, at least for the above choice of scales. Only for negative ?7iab where the cross section 
changes sign are the corrections large. 

An important reason why inclusion of the NLO corrections is generally vital is that they 
are expected to considerably reduce the dependence of the cross sections on the unphysical 
factorization and renormalization scales. This improvement in scale dependence when going 
from LO to NLO is a measure of the impact of the NLO corrections, the residual scale depen- 
dence at NLO perhaps providing a rough estimate of the relevance of even higher order QCD 
corrections. Figure El shows the scale dependence of the spin-dependent cross section at LO 
and NLO. In each case the shaded bands indicate the uncertainties from varying the unphysical 
scales in the range Pt/2 < /^r = /^/ = A*/ < '^Vt- The solid lines are for the choice where all 
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Figure 4: NLO direct and resolved photon contributions in the MS and DIS^ factorization schemes. 
Also shown is their sum which has to be independent of theoretical conventions. 



scales are set to p^. One can see that the scale dependence indeed becomes much smaller at 
NLO and that the lowest order approximation may capture the main features of the process 
but will not provide a satisfactory quantitative understanding. 

Figure m demonstrates, as discussed in Section II, that the direct and resolved contributions 
individually depend on the choice of the factorization scheme, but that the physical, i.e., ex- 
perimentally accessible, cross section does not. To make this point, we transform the polarized 
cross sections from the MS to the DIS^ scheme as described in Eqs. (pT]|) - (fT^ . 

A first determination of the parton content of circularly polarized photons could be one of 
the key physics goals of the spin physics program at a future lepton-proton collider. Figure El 
shows the experimentally relevant spin asymmetry 



da dijA 



daA 



(15) 



for tt'^ photoproduction at eRHIC. As before we have integrated over pt > 4 GeV and have 
chosen all scales to be pt- At positive rapidities the spin asymmetry shows the expected 
dependence on the choice of the polarized photon densities. In order to be able to judge whether 
a future measurement can resolve such a difference we also give in Figure El an estimate for the 
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Figure 5: Spin asymmetry for 7r° photoproduction in NLO QCD for the two extreme sets of polarized 
photon densities and two different choices of spin-dependent proton distributions. The "error bars" 
indicate the statistical accuracy anticipated for eRHIC assuming an integrated luminosity of lfb~^. 



expected statistical errors: 



i] 



(16) 



where Ve and Vp are the polarization of the lepton and proton beam, respectively, C denotes 
the integrated luminosity of the collisions, and o"bin is the unpolarized cross section integrated 
over bins in r^iab- We have used Ve,p = 0.7 and £ = 1/fb, which are targets for eRHIC. It should 
be noted that an integrated luminosity of £ = 1/fb is expected to be accumulated after only a 
few weeks of running so that the statistical accuracy may eventually be much better than the 
one shown in Fig. |S1 

To demonstrate that the sensitivity to A/"^ is not masked by the dependence on the spin- 
dependent proton densities, we show in Fig. |S] also the results obtained for the GRSV "max. 
Agf" set [16j, which is characterized by a much larger gluon density A (7 than the one in the 
GRSV "standard" set. It is evident that at large positive rapidities the asymmetry is indeed 
mainly determined by the polarized photon structure. In the region of ?7iab < a measurement 
of would be also an excellent source of information on the polarized gluon density in the 
nucleon. We note that, in any case, by the time eRHIC is commissioned we expect the spin 
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Figure 6: Relative contributions of different partonic scatterings ab cX to the NLO (MS) 
resolved photon cross section for 7r° production at eRHIC, using the "standard" set of GRSV [16j 
for the proton and the "maximal" scenario for the spin-dependent photon densities. 

structure of the nucleon to be known with much better accuracy than at present. 

Figure IHl illustrates how different partonic scatterings ab cX contribute to the NLO 
(MS) polarized resolved cross section, versus rapidity. Recall that parton "a" emerges from 
the photon and parton "6" from the proton. Qualitatively similar results are obtained for the 
"minimal" photon scenario and hence not shown here. 

IV Conclusions 

In this paper we have presented the first complete NLO calculation for the spin-dependent 
photopro duct ion of inclusive high-p^^ hadrons, consistently including the direct and resolved 
photon contributions at NLO. Our results are relevant for future polarized lepton-proton collider 
projects like eRHIC, to which we have tailored our numerical results. The spin asymmetry for 
this process is a promising tool for a first determination of the spin structure of circularly 
polarized photons. We found that the QCD corrections to the polarized cross section are under 
good control. The polarized cross section shows a significant reduction of scale dependence 
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when going from LO to NLO. 



Wc finally note that a single process alone, like pion photoproduction, will allow to establish 
the very existence of a polarized resolved-photon contribution and may perhaps point towards 
one of the two scenarios discussed above. It will, however, not be sufficient to fully disentangle 
the various quark flavor distributions and the gluon density of the polarized photon. This 
will require, as for the nucleon (spin) structure, a multitude of processes, combined in a global 
QCD analysis. Of course, at eRHIC other photoproduction processes, like jets, di-jets, or heavy 
flavors, will add further valuable information. Measurements at a future polarized linear e'^e~ 
collider would be also extremely helpful in this respect. 
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